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Introduction 



In the Standard Model (SM) the decays of pseudoscalar mesons to light leptons are helicity 
suppressed. In particular, the SM width of P^ —> l^v decays with P = it, K, D, B (denoted P& 
in the following) is 

where Gf is the Fermi constant, Mp and Mi are meson and lepton masses, fp is the decay 
constant, and V qq , is the corresponding Cabibbo-Kobayashi-Maskawa matrix element. Although 
the SM predictions for the Pi2 decay rates are limited by hadronic uncertainties, their specific 
ratios do not depend on fp and can be computed very precisely. In particular, the SM pre- 
diction for the ratio Rk = T(K e 2) /T^K^) of kaon leptonic decay widths inclusive of internal 
bremsstrahlung (IB) radiation is pQ 

= (0 (|N|) 2 < 1 + W ^ = P™^) >< P) 

where SRqed = (— 3.79 ± 0.04)% is an electromagnetic correction due to the IB and structure- 
dependent effects. 

Within certain two Higgs doublet models (2HDM), including the minimal supersymmetric 
model (MSSM), Rk is sensitive to lepton flavour violating (LFV) effects appearing at the one- 
loop level via the charged Higgs boson (H^) exchange [21 [3], representing a unique probe into 
mixing in the right-handed slepton sector [JJ. The dominant contribution due to the LFV 
coupling of the is 



nLFV pSM 



1+ rifL lA^Ptan^ 



M K y ( M T \ 2 IA 3i |2 , 



(3) 

where tan f3 is the ratio of the two Higgs vacuum expectation values, and |A^| is the mixing 
parameter between the superpartners of the right-handed leptons, which can reach ~ 10 -3 . 
This can enhance Rk by 0(1%) without contradicting any experimental constraints known at 
present, including upper bounds on the LFV decays r — > eX with X = rj, 7, fxfi. On the other 
hand, Rk is sensitive to the neutrino mixing parameters within the SM extension involving a 
fourth generation [5]. 

The first measurements of Rk were performed in the 1970s [6l O [8] ; the current PDG world 
average [9] is based on a more precise recent result [10] Rk = (2.493 ± 0.031) x 10~ 5 . A new 
measurement of Rk based on a part of the data sample collected by the NA62 experiment at 
CERN in 2007 is reported in this letter. The analyzed K e 2 sample is ~ 4 times larger than the 
total world sample, allowing a measurement of Rk with a precision well below 1%. 



1 Beam, detector and data taking 

The beam line and setup of the NA48/2 experiment I12j have been used for the NA62 data 
taking in 2007. Experimental conditions have been optimized for the K^jK^ measurement. 

1.1 Kaon beam 

The beam line was originally designed to deliver simultaneous unseparated K + and K~ beams 
derived from the primary 400 GeV/c protons extracted from the CERN SPS. In 2007, the muon 
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sweeping system was optimized for the positive beam, and the sample used for the present 
analysis was collected with the K + beam only. Positively charged particles within a narrow 
momentum band with a central momentum of 74.0 GeV/c and a spread of 1.4 GeV/c (rms) 
are selected by the first two magnets in a four dipole achromat and by momentum-defining 
slits incorporated into a 3.2 m thick copper /iron proton beam dump, which also blocks the 
negatively-charged particles. The beam subsequently passes through acceptance-defining and 
cleaning collimators and a set of four quadrupoles of alternating polarity, as well as muon 
sweeping magnets, before entering the fiducial decay volume contained in a 114 m long cylindrical 
vacuum tank with a diameter of 1.92 m upstream, increasing to 2.4 m downstream. 

With about 1.8 x 10 12 primary protons incident on the target per SPS pulse of about 4.8 s 
duration repeating every 14.4 or 16.8 s, the secondary beam flux at the entrance of the decay 
volume is 2.5 x 10 particles per pulse. The fractions of K + , tt + , p + , e + and /i + in the secondary 
beam are 0.05, 0.63, 0.21, 0.10 and 0.01, respectively. The fraction of those beam kaons decaying 
in the vacuum tank at nominal momentum is 18%. The beam transverse size at the entrance to 
the decay volume is 5x = 5y = 4 mm (rms), and its horizontal and vertical angular divergences 
are about 20 /xrad (rms). The beam central momentum, transverse position at the entrance to 
the vacuum tank and direction varied slowly over time with respect to the nominal ones in the 
ranges of ~ 0.1 GeV/c 2 , ~ 1 mm and ~ 10 /^rad, respectively. 

The beam line also transmits certain off-momentum charged kaons and pions punching 
through the proton beam dump with a suppression factor of ~ 10~ 3 . However the subsequent 
decays of these particles do not contribute to the signal region of the present analysis. 

1.2 Detector 

The charged particle properties are measured in a magnetic spectrometer, housed in a tank 
filled with helium at nearly atmospheric pressure, placed after the decay volume and separated 
from the vacuum by a thin (~ 0.4% radiation lengths Xq) Kevlar® window. The spectrometer 
comprises four drift chambers (DCHs) [13], two upstream and two downstream of a dipole 
magnet which gives a horizontal transverse momentum kick of 265 MeV/c to singly-charged 
particles. Each DCH is composed of eight planes of sense wires, and provides a spatial resolution 
of 90 /um in each projection. The measured momentum resolution is cr p /p = 0.48% © 0.009% • p, 
where p is expressed in GeV/c. 

A plastic scintillator hodoscope (HOD) producing fast trigger signals and providing precise 
time measurements of charged particles is placed after the spectrometer. It consists of a plane 
of vertical strips, followed by a similar plane of horizontal strips (128 counters in total). Both 
planes have regular octagonal shapes and a central hole for the passage of the beam. 

The HOD is followed by a quasi-homogeneous liquid krypton electromagnetic calorimeter 
(LKr) [T5] used for lepton identification and as a photon veto in the present analysis. The LKr 
is 127 cm (or 27Xo) thick along the beam, with projective readout consisting of copper /beryllium 
ribbons extending from the front to the back of the detector. The 13248 readout cells have a 
transverse size of approximately 2x2 cm 2 each and have no longitudinal segmentation. The 
energy resolution is a E /E = 0.032/\/£ © 0.09/E © 0.0042 (E in GeV). The spatial resolution 
for the transverse coordinates x and y of an isolated electromagnetic shower is a x = a y = 
0.42/ 'VE © 0.06 cm (E in GeV). 

An aluminium beam pipe of 158 mm outer diameter and 1.1 mm thickness traversing the 
centres of all detector elements allows the undecayed beam particles to continue their path in 
vacuum. The outer transverse sizes of the subdetectors are about 2.4 m. 
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1.3 Trigger and data acquisition 

A minimum bias trigger configuration has been employed, resulting in high efficiency. The K e2 
trigger condition consists of coincidences of signals in the two HOD planes (the Q\ signal), loose 
lower and upper limits on DCH hit multiplicity (the 1-track signal), and LKr energy deposit 
(-E-LKr) of at least 10 GeV. The K^ 2 trigger condition requires a coincidence of the Q± and 1-track 
signals downscaled by a factor D = 150. The non-downscaled K^ 2 trigger rate is 0.5 MHz, and 
is dominated by beam halo muons; the K e2 trigger rate is about 10 kHz. Downscaled control 
samples based on trigger signals from the DCHs, HOD and LKr have been collected to monitor 
the performance of the main trigger signals. The data taking took place during four months 
starting in June 2007. About 40% of the 350k recorded good SPS spills are used for the present 
analysis. 

2 Analysis strategy 

The analysis strategy is based on counting the numbers of reconstructed K e2 and K^ 2 candidates 
collected concurrently. Therefore the analysis does not rely on an absolute beam flux measure- 
ment, and several systematic effects (due to beam simulation, accidental activity, charged track 
reconstruction, Q\ trigger efficiency, and time-dependent effects) cancel at first order. 

Due to the significant dependence of acceptance and background on lepton momentum, the 
Rk measurement is performed independently in 10 momentum bins covering a range from 13 to 
65 GeV/c. The lowest momentum bin spans 7 GeV/c, while the others are 5 GeV/c wide. The 
selection criteria have been optimized separately in each momentum bin. The data samples in 
the momentum bins are statistically independent, however the systematic errors are partially 
correlated. The ratio Rk in each bin is computed as 

1_ N(K e2 )-N B (K e2 ) A{K^ 2 ) fji x e(K^) 1 
K D ' N(K^ 2 ) - iV B (i^ 2 ) ' A(K e2 ) ' f e x e(K e2 ) ' / LKr ' 1 ' 

where N{K^ 2 ) are the numbers of selected Kg} candidates (£ = e,fx), N-q{K^ 2 ) are the numbers 
of background events, A{K^ 2 ) / A(K e2 ) is the ratio of the geometric acceptances (referred to as 
the acceptance correction in the following), fg are the lepton identification efficiencies, e(Ki 2 ) 
are the trigger efficiencies, /LKr is the global efficiency of the LKr readout (which affects only 
the K e2 selection), and D = 150 is the K^ 2 trigger downscaling factor. 

To evaluate the acceptance correction and the geometric parts of the acceptances for back- 
ground processes entering the computation of Nb(Ki 2 ), a detailed Monte Carlo (MC) simulation 
based on Geant3 [15] is used. It includes a description, with time variations, of the beam line 
optics, the full detector geometry, materials, magnetic fields, local inefficiencies of DCH wires, 
and inactive LKr cells (0.8% of channels). Particle identification, trigger and readout efficiencies 
are measured directly from data. 

3 Data analysis 

3.1 Event reconstruction and selection 

Charged particle trajectories are reconstructed from hits and drift times in the spectrometer. 
Track momenta are evaluated using a detailed magnetic field map. Fine calibrations of spec- 
trometer field integral and DCH alignment are performed by monitoring the mean reconstructed 
K + 7r + 7r + 7r~ invariant mass, and the missing mass in K^ 2 decays. 

Clusters of energy deposition in the LKr are found by locating the maxima in the digitized 
pulses from individual cells in both space and time and accumulating the energy within a radius 
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of approximately 11 cm. Shower energies are corrected for energy outside the cluster boundary, 
energy lost in inactive cells and cluster energy sharing. The energy response has been calibrated 
with samples of positrons from K + — > decays. 

Due to the topological similarity of K e 2 and decays, a large part of the selection is com- 
mon for the two decay modes, which leads to significant cancellations of the related systematic 
uncertainties. The main selection criteria are listed below. 

• Exactly one reconstructed charged particle track (lepton candidate) geometrically consis- 
tent with originating from a kaon decay is required. The electric charge of the track must 
be positive. 

• The extrapolated track impact points in the DCHs, HOD and LKr must be within their 
geometrical acceptances. The LKr acceptance condition includes appropriate separations 
from the detector edges and inactive cells. 

• The reconstructed track momentum must be in the range 13 to 65 GeV/c. The lower limit 
ensures high efficiency of the -Elkt > 10 GeV trigger condition. Above the upper limit, 
the analysis is affected by large uncertainties due to background subtraction. 

• No LKr clusters with energy E > E Yeto = 2 GeV and in time with the track are allowed, 
unless they can be associated to the track via direct energy deposition or bremsstrahlung. 
(Most clusters due to bremsstrahlung in front of the spectrometer magnet are resolved 
from those directly deposited by the track). This requirement provides a photon veto for 
suppression of backgrounds from K + — > e + uj, K + — > 7r°e + v, and K + — > 7r + 7r° decays. 
However the veto is not hermetic due to the beam pipe and the limited transverse size of 
the LKr. 

• The decay vertex is reconstructed as the point of closest approach of the lepton candidate 
track extrapolated upstream, and the kaon beam axis. The measured stray magnetic field 
in the vacuum tank is taken into account. The position of the kaon beam axis is monitored 
with a sample of fully reconstructed K + — > tt + tt + tt~ decays. 

• The distance from the kaon decay vertex to the beginning of the vacuum tank is re- 
quired to exceed a minimum value ranging from 8 m at low lepton momentum to 43 m 
at high momentum, which removes the bulk of the beam halo background (discussed in 
Section 

• For further suppression of the beam halo and several other backgrounds, the reconstructed 
closest distance of approach of the track to the beam axis must not exceed 3.5 cm. 

The following two main criteria are used to distinguish from decays. 

• The kinematic identification of K e 2 (K^) decays is based on constraining the reconstructed 
squared missing mass in the positron (muon) hypothesis: 

- Ml < Ml iss (£) = (P K - P t f < Ml (5) 

Here Pk is the average kaon four-momentum (monitored in time with K + — > tt + tt + tt^ 
decays), and Pi is the reconstructed lepton four-momentum (under the positron or muon 
mass hypothesis). The limits M 2 and M 2 have been optimized for each lepton momentum 
bin, taking into account the M^ iss (£) resolution (which varies from 0.0025 (GeV/c 2 ) 2 at 
mid track momentum to 0.005 (GeV/c 2 ) 2 at low and high track momentum), the radiative 
mass tails, and the background conditions. M 2 varies between 0.013 and 0.016 (GeV/c 2 ) 2 
and Mf between 0.010 and 0.013 (GeV/c 2 ) 2 . The kinematic separation of K e 2 and 
decays is illustrated in Fig. [T£i. 
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Figure 1: (a) Squared missing mass assuming the positron mass hypothesis M^ iss (e) as a function 
of lepton momentum for reconstructed K e 2 and decays (data); (b) K e 2 vs kinematic 
separation (standard deviations) as a function of lepton momentum; (c) E/p spectra of positrons 
and muons (data); the positron identification limits for p > 25 GeV/c are indicated by arrows. 



• The lepton identification is based on the ratio E/p of energy deposition in the LKr to 
momentum measured by the spectrometer. Charged particles with (E/p) m \ n < E/p < 1.1, 
where {E/p) m \ n = 0.95 for p > 25 GeV/c and (E/p) m \ n = 0.9 otherwise, are identified 
as positrons. At low lepton momenta, the background from particle mis-identification is 
negligible. For p > 25 GeV/c, the larger (E/p) m [ n limit minimises the net uncertainty from 
K^2 background subtraction and particle mis-identification inefficiency. Charged particles 
with E/p < 0.85 are classified as muons. The data E/p spectra of positrons and muons 
are shown in Fig. [lb. 

3.2 The K e2 sample 

The number of K e 2 candidates in the signal region is N(K e 2) = 59813. The sources of back- 
ground in the K e 2 sample are discussed below. 

3.2.1 background 

Kinematic separation of K e 2 from K^i decays is achievable at low lepton momentum only (p < 
35 GeV/c), as shown in Figs. [Tb, and [lb. At high lepton momentum, the decay with a 
mis-identified muon (E/p > 0.95, see Fig. [lb) is the largest background source. The dominant 
process leading to mis-identification of the muon as a positron is 'catastrophic' bremsstrahlung 
in or in front of the LKr leading to significant energy deposit in the LKr. Mis-identification 
due to accidental LKr clusters associated with the muon track is negligible, as concluded from 
a study of the sidebands of track-cluster time difference and distance distributions. 

The muon mis-identification probability P^ e has been measured as a function of momentum. 
To collect a muon sample free from the typical ~ 10 -4 positron contamination due to fj, — > e 
decays, a 9.2Ao thick lead (Pb) wall covering ~ 20% of the geometric acceptance was installed 
approximately 1.2 m in front of the LKr calorimeter (between the two HOD planes) during a 
dedicated period of data taking with K + and K~ beams. The K e 2 sample collected with the 
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Figure 2: (a) Mis-identification probability for muons traversing the lead wall, Pnei f° r 
(E/p) m i n = 0.95 as a function of momentum: measurement (solid circles with error bars) 
and simulation (solid line), (b) Correction factors /pb = P^e/P^e f° r the considered values 
of (E/p) m [ Q , as evaluated with simulation. Dotted lines in both plots indicate the estimated 
systematic uncertainties of the simulation. 

Pb wall installed is not used for the Rk measurement. The component from positrons which 
traverse the Pb wall and are mis-identified as muons from decay with p > 30 GeV/c and 
E/p > 0.95 is suppressed down to a negligible level (~ 10~ 8 ) by energy losses in the Pb. 

However, muon passage through the Pb wall affects the measured Pj^ via two principal 
effects: 1) ionization energy loss in Pb decreases and dominates at low momentum; 2) 
bremsstrahlung in Pb increases P^ e and dominates at high momentum. To evaluate the correc- 
tion factor /pb = P/_ ie /P^e, a dedicated MC simulation based on Geant4 (version 9.2) [16] has 
been developed to describe the propagation of muons downstream from the last DCH, involving 
all electromagnetic processes including muon bremsstrahlung |17j . 

The measurements of PY~ in momentum bins compared with the results of the MC simulation 
and the correction factors /pb obtained from simulation, along with the estimated systematic 
uncertainties of the simulated values, are shown in Fig. [2j The relative systematic uncertainties 
on P^e and Pj^ obtained by simulation have been estimated to be 10%, and are mainly due to the 
simulation of cluster reconstruction and energy calibration. However the error of the ratio /pb = 
P^e/Pjie i s significantly smaller (<5/pb//pb = 2%) due to cancellation of the main systematic 
effects. The measured is in agreement with the simulation within their uncertainties. 

The positive correlation between the reconstructed M^j (e) and E/p, which are both com- 
puted using the reconstructed track momentum, leads to an apparent dependence of P^ e on 
M^ iss (e). This effect is significant for intermediate lepton momenta where the background 
comes from events with underestimated M^ iss (e) and a smaller muon mis-identification proba- 
bility (see Fig. QJi). This correlation has been taken into account. 

The K^2 background contamination integrated over lepton momentum has been computed 
to be (6.11 ± 0.22)% using the measured Pj^ corrected by /pb- The quoted error comes from 
the limited size of the data sample used to measure P^ e b (0.16%), the uncertainty <5/pb (0.12%), 
and the model-dependence of the correction for the M^ iss (e) vs E/p correlation (0.08%). The 
first error component is uncorrelated between the lepton momentum bins, while the others are 
fully correlated. 
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As a stability check, the evaluation of P^ e has been performed with an additional requirement 
that the energy deposit in the HOD counters downstream from the Pb wall is small (limited 
to the equivalent of 1.5 to 3 minimum ionizing particles), which strongly suppresses muons 
undergoing bremsstrahlung in the Pb wall. The stability of P^ e is consistent with the assigned 
uncertainty <5/pb- Additionally, a stability check of Rk with respect to variation of (E/p) m in 
for lepton momentum p > 25 GeV/c in the range from 0.90 to 0.97 has been performed. The 
observed relative stability of Rk within ±0.2%, although the background varies from 17% 
to 3%, is consistent with the uncertainty assigned to the background. 

The K^2 decay also contributes to background via /x + — > e + vv decays in flight. Energetic 
forward secondary positrons compatible with K e 2 kinematics and topology are suppressed by 
muon polarisation effects [18 . Radiative corrections to the muon decay [19] lead to a further 
~ 10% relative background suppression. This background contamination has been estimated to 
be (0.27 ± 0.04)%, where the dominant uncertainty is due to the simulated statistics. 

3.2.2 K + — > background 

Rk is defined to be fully inclusive of internal bremsstrahlung (IB) radiation pQ. The structure- 
dependent (SD) K + — > e + z^7 process [20], [21] may lead to a K e 2 signature if the positron is 
energetic and the photon is undetected. In particular, the SD + component with positive photon 
helicity peaks at high positron momentum in the K + rest frame {E* ~ M^-/2) and has a similar 
branching ratio to K e 2- The background due to K + — > e + u^f (SD - ) decay with negative photon 
helicity peaking at E* ~ and the interference between the IB and SD processes are 

negligible. 

The SD + background contribution has been estimated by MC simulation as (1.07 ± 0.05)%, 
using a recent measurement of the K + — > e + v~/ (SD + ) differential decay rate [10j . The quoted 
uncertainty is due to the limited precision on the form factors and decay rate, and is therefore 
correlated between lepton momentum bins. A stability check of Rk against variation of the 
E vcto limit in a wide range has been performed. While the K + — > e + uj (SD + ) background is 
enhanced by a factor of 4.5 for E veto = 14 GeV with respect to E vcto = 2 GeV, Rk remains 
stable within ±0.1%, which is consistent with the above uncertainty. 

3.2.3 K+ 7v°e+u and K~ — >- tv+tt backgrounds 

The K + —7- 7r°e + i/ decay produces a K e 2 signature if the only reconstructed particle is an e + 
from K + or tt° Dalitz (vr^, — > 7e + e~) decays. The K + — > 7r + 7r° decay leads to a K e 2 signature 
if the only reconstructed particle is a ir + mis-identified as e + , or an e + from a — > je + e~ 
decay. The pion mis-identification probability (0.95 < E/p < 1.1) has been measured to be 
(0.41 ± 0.02)% in the relevant momentum range from samples of K + — > tt + tt° and — > -K^e^v 
decays (the latter collected during a special run). 

Kinematically, K + — > -K®e + v and K + — > 7r + 7r° decays can be reconstructed with low miss- 
ing mass in the K e 2 signal region, either because the charged track undergoes a large multiple 
scatter or because the kaon is in the high momentum tail of the beam distribution. The system- 
atic uncertainties due to subtraction of these backgrounds have been estimated as 50% of the 
contributions themselves, due to the limited precision of the simulation of the kaon momentum- 
distribution tails. The backgrounds are at a level below 0.1%. 

3.2.4 Beam halo background 

As no tracking is available in the beam region to tag an incoming kaon, beam halo muons can 
become a source of background to K e 2 decays in case of /x + — > e + v e r , ll decay or muon mis- 
identification as a positron. The choice of the signal region in terms of the longitudinal position 
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Table 1: Summary of backgrounds in the K e 2 sample. 

Source N B /N(K e2 ) 

K^ 2 (6.11 ± 0.22)% 

K^^e) (0.27 ±0.04)% 

K+ -»• e+v-f (SD+) (1.07 ± 0.05)% 

K+ -> 7r°e+z/ (0.05 ± 0.03)% 

K + -> vr+vr (0.05 ± 0.03)% 

Beam halo (1.16 ±0.06)% 

Total (8.71 ± 0.24)% 



of the kaon decay vertex has been dictated by the kinematic distribution of this background 
(which peaks in the upstream part of the vacuum volume). 

The halo background has been measured directly by reconstructing the K^ 2 candidates from 
one control data sample collected with the K~ beam transmitted by the beam line and the 
K + beam (but not its halo) blocked, and another control data sample collected with both 
K + and K~ beams blocked. The control samples are normalised to the data in the region 
—0.3 < M^ liss ( / u) < —0.1 (GeV/c 2 ) 2 populated mainly by beam halo events. The 'cross-talk' 
probability to reconstruct a K^ 2 candidate due to a K~ decay with e + emission (K~ — > -k® d I~v, 
K~ — > ^iTp, K~ —7- £~ve + e~ with i = e,fi) is at the level of ~ 10 -4 and is taken into 
account. The halo background rate and kinematical distribution are qualitatively reproduced 
by a simulation of the beam line. 

The background contamination has been estimated to be (1.16 ± 0.06)%, where the error 
comes from the limited size of the control samples (uncorrelated between lepton momentum 
bins) and the normalisation uncertainty due to decays of beam kaons and pions upstream of the 
decay volume (correlated between momentum bins). 

3.2.5 Summary of backgrounds in the K e 2 sample 

Backgrounds integrated over lepton momentum are summarised in Table [TJ The total back- 
ground contamination is (8.71 ± 0.24)%; its uncertainty is smaller than the relative statistical 
uncertainty of 0.43%. The M^ iss (e) and lepton momentum distributions of K e 2 candidates and 
backgrounds are shown in Fig. [3j 

3.3 The Kfj_2 sample 

The number of K^ 2 candidates collected with a trigger chain involving downscaling by a factor 
of 150 is N{K^2) — 1-803 x 10 7 . The only significant background source in the sample is 
the beam halo. Its contribution is mainly at low muon momentum, and has been measured to 
be (0.38 ± 0.01)% using the same technique as for the K e2 sample. The M^ aiss (fi) and muon 
momentum spectra of K^ 2 candidates and the halo background are presented in Fig. HI 

3.4 Geometrical acceptance correction 

The ratio of geometric acceptances A(K^ 2 )/A(K e 2) in each lepton momentum bin has been 
evaluated with MC simulation. The radiative K + — > e + i^7 (IB) process, which is responsible for 
the loss of about 5% of the K e 2 acceptance by increasing the reconstructed M^ iss (e), is taken 
into account following [20] , with higher order corrections according to [22l [23] . 
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Figure 3: (a) Reconstructed squared missing mass M^ iss (e) distribution of the K e 2 candidates 
compared with the sum of normalised estimated signal and background components. The small 
discrepancy between data and MC at low M^ iss (e) is due to the limited precision of MC beam 
description, which is taken into account by systematic uncertainty due to the acceptance correc- 
tion, (b) Lepton momentum distributions of the K e % candidates and the dominant backgrounds; 
the backgrounds are scaled for visibility. 




M^ iss (n), (GeV/c 2 ) 2 Lepton momentum, GeV/c 



Figure 4: (a) Reconstructed squared missing mass M^ iss (/i) distribution of the K^2 candidates 
compared with the sum of normalised estimated signal and background components. The deficit 
of reconstructed MC events in the region of the K + — > tt + tt° peak is due to the limited pre- 
cision of the beam simulation, and is mostly outside the signal region, (b) Lepton momentum 
distributions of the candidates and the beam halo background (the latter is scaled for 
visibility) . 
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Lepton momentum, GeV/c Positron momentum, GeV/c 

Figure 5: (a) The acceptance correction A(Kf Jj 2)/A(K e 2) in lepton momentum bins; the correc- 
tions neglecting internal (IB) and external (EB) bremsstrahlung radiation are also presented, 
(b) The measured positron identification inefficiency 1 — f e in lepton momentum bins; uncer- 
tainties in bins are partially correlated. The lower inefficiency for p < 25 GeV/c is due to the 
relaxed positron identification requirement discussed in Section 13.11 

The acceptance correction is strongly influenced by bremsstrahlung suffered by the positron 
in the material upstream of the spectrometer magnet (Kevlar window, helium, DCHs). This 
results in an almost momentum-independent loss of K e 2 acceptance of about 6%, mainly by 
increasing the reconstructed M^ iss (e). The relevant material thickness has been measured by 
studying the spectra and rates of bremsstrahlung photons produced by low intensity 25 GeV/c 
and 40 GeV/c electron and positron beams steered into the DCH acceptance, using special data 
samples collected in the same setup by the NA48/2 experiment in 2004 and 2006. Using these 
measurements, the material thickness during the 2007 run has been estimated to be (1.56 ± 
0.03)%A"o. The quoted uncertainty is dominated by the limited knowledge of helium purity in 
the spectrometer tank; its measured purity of (92 ± 4)% corresponds to a thickness of (0.26 ± 
0.03)%Ao. This translates into a systematic uncertainty on Rk- 

The acceptance correction A(K ^2) / 'A(K e 2) in lepton momentum bins is presented in Fig. [5k.. 
The corrections evaluated without internal (IB) and external (EB) bremsstrahlung radiation 
are also presented to illustrate the magnitudes of the corresponding effects. The correction is 
enhanced at low lepton momentum because the radial distributions of positrons from K e 2 decays 
in the DCH planes are wider than those of muons from decays, and low momentum leptons 
are not fully contained within the geometric acceptance due to the limited transverse sizes of 
the DCHs. 

The track reconstruction inefficiency due to interactions in spectrometer material is included 
into the acceptance correction. Simulation of the positron track reconstruction inefficiency 
(which is ~ 10 -3 in the analysis track momentum range) has been validated with a sample of 
K+ -> 7r + 7r2) decays. The muon track reconstruction inefficiency evaluated with MC simulation 
is ~ 2 x 10~ 4 . Systematic effects due to imperfect simulation of the reconstruction efficiency are 
negligible. 

Apart from helium purity, the main sources of systematic uncertainty of the acceptance 
correction are the limited knowledge of beam profile and divergence, and the simulation of soft 
radiative photons. A separate uncertainty has been assigned to account for the finite precision 
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of the DCH alignment. 



3.5 Lepton identification efficiencies 

The E/p ratio provides powerful particle identification criteria. The momentum-dependent 
positron identification window (E/p) m [ n < E/p < 1.1 includes more than 99% of the K e 2 
events, while suppressing muons by a factor of 1/P^ e ~ 10 6 . The requirement E/p < 0.85 leads 
to a negligible inefficiency of the muon identification. 

A pure sample of 4 x 10 7 positrons, selected kinematically from K + — > 7r°e + z/ (charged 
K e z) decays collected with the K e 2 trigger concurrently with the main fQ 2 data set, is used to 
calibrate the energy response of each LKr cell and to study f e with respect to local position 
and time stability. However, the momentum range of the positrons from charged K e 3 decays 
is kinematically limited, preventing a sufficiently precise measurement of f e above 50 GeV/c. 
Therefore a dedicated data sample was recorded in a special 15 hour long run with a broad 
momentum band K® beam. Electrons and positrons from the 4 x 10 6 collected K® — > 7r ± e =F ^ 
(neutral K e %) decays allow the determination of f e in the whole analysis momentum range. 

The measurements of f e have been performed in bins of lepton momentum; a finer binning is 
used inside the lowest of the 10 standard bins to improve the determination of local inefficiencies, 
which peak at low momentum. Separate measurements have been performed for several identified 
groups of LKr cells with higher local inefficiencies. Efficiency measurements with the charged 
and neutral K e ^ decays agree to better than 0.1%. Fig. [5b shows the measurements of 1 — 
f e in momentum bins used to evaluate corrections to Rk, obtained as the weighted mean of 
charged and neutral K e ^ measurements for momenta up to 50 GeV/c, and as neutral kaon 
measurements for higher momenta. The inefficiency averaged over the K e 2 sample is 1 — f e = 
(0.73 ± 0.05)%, where the uncertainty takes into account the statistical precision and the small 
differences between charged and neutral kaon results. 

3.6 Trigger and readout efficiencies 

The efficiency of the Q\ trigger condition has been measured using events triggered with 
a control LKr signal. The inefficiency integrated over the Ku2 sample is (1.4 ± 0.1)%. As a 
consequence of its geometric uniformity and the similarity of the K e 2 and distributions over 
the HOD plane, it nearly cancels between the K e 2 and samples, and the residual systematic 
bias is negligible. The inefficiency of the 1-track condition also largely cancels in the ratio Rk, 
but is anyway negligible. 

Thus the trigger efficiency correction e{K^2) / e {Ke2) is determined by the efficiency e(-ELKr) 
of the LKr energy deposit trigger signal -Elkj > 10 GeV entering the K e 2 trigger chain only. The 
inefficiency 1 — e(i?LKr) is only significant in the lowest lepton momentum bin of (13, 20) GeV/c, 
which is close to the trigger energy threshold and is thus affected by the online energy resolution. 
A sample of events triggered with a control Q\ signal passing all K e 2 selection criteria except the 
M^ iss (e) constraint, therefore dominated by K e ^ events with two lost photons, has been used to 
measure 1 — e(£ , LKr) in the lowest momentum bin to be (0.41 ± 0.05 s t a t.)%- Corrected for the 
difference of positron distributions in the LKr plane between the K e 2 sample and the control 
sample, it translates into 1 — e(i?LKr) = (0.61 ± 0.20)% for the K e 2 sample. The correction 
and its uncertainty are significant due to the presence of several locally inefficient regions. The 
resulting uncertainty on Rk is negligible. 

Energetic photons not reconstructed in the LKr may initiate showers by interacting in the 
DCHs or the beam pipe, causing the DCH hit multiplicities to exceed the limits allowed by 
the 1-track trigger condition. Among the backgrounds, only the K + — > e + uj (SD + ) receives a 
non-negligible correction due to the 1-track inefficiency. The inefficiency for K + — > 7r°e + f events 
with two lost photons has been measured to vary in the range from 0.1 to 0.3 depending on 
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Figure 6: Measurements of Rk in lepton momentum bins with their uncorrelated statistical 
uncertainties and the partially correlated total uncertainties. The average Rk and its total 
uncertainty are indicated by a band. 



track momentum. The extrapolation of this result to K + — > e + w~y (SD + ) with one lost photon 
relies on simulation. The corresponding uncertainty has been propagated into Rk. 

The global LKr readout inefficiency, affecting the K e 2 reconstruction only, has been measured 
using an independent readout system to be 1 — /LKr = (0.20 ± 0.03)%, stable in time. This 
measurement has been cross-checked, with limited precision, by a study of the LKr response in 
a sample of tt^ d — > Ae^ 1 decays reconstructed from spectrometer information only. 

4 Result and discussion 

A x 2 fit t° the measurements of Rk in the 10 lepton momentum bins has been performed, taking 
into account the bin-to-bin correlations between the systematic errors. To validate the assigned 
systematic uncertainties, extensive stability checks have been performed in bins of kinematic 
variables and by varying selection criteria and analysis procedures. The fit result is 

R K = (2.487 ± 0.01W ± 0.007 syst .) x 10" 5 = (2.487 ± 0.013) x 10 -5 , (6) 

with x 2 /ndf = 3.6/9. The individual measurements with their statistical and total uncertainties, 
and the combined result are displayed in Fig. [6J The uncertainties of the combined result are 
summarised in Table [2j 

This is the most precise Rk measurement to date. It is consistent with the KLOE measure- 
ment [10] and the SM expectation R^ = (2.477 ± 0.001) x 10 -5 , and can be used to constrain 
multi-Higgs [2] and fourth generation [5] new physics scenarios. The experimental accuracy is 
still an order of magnitude behind the SM accuracy, which motivates further precision measure- 
ments of Rk- 
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Table 2: Summary of the uncertainties on Rk- 



Source 


SR K x 10 5 


Statistical 


0.011 


background 


0.005 


K + —7- e + uj (SD + ) background 


0.001 


K + — > ir e + v, K + — > 7r + 7r° backgrounds 


0.001 


Beam halo background 


0.001 


Helium purity 


0.003 


Acceptance correction 


0.002 


Spectrometer alignment 


0.001 


Positron identification efficiency 


0.001 


1-track trigger efficiency 


0.002 


LKr readout inefficiency 


0.001 


Total systematic 


0.007 


Total 


0.013 
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